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PREFACE 


The subject “Recursion Sequences” is related with the secon- 
dary school curriculum (arithmetic and geometric progres- 
sions, the sequence of squares of natural numbers, the sequence 
of coefficients in the quotient of two polynomials in 
ascending powers etc.). It is at the same time a real 
though a small mathematical theory*, complete and clear 
like everything created by the eminent masters of analysis 
who founded it. 

The fundamentals of the theory of recursion sequences 
were worked out and published in the twenties of the 
XVIII-th century by the French mathematician Abraham 
De Moivre [the formula (cosa-+isina)"=cosna-+isinna 
bears his name] and by one of the first members of the 
Petersburg Academy of Sciences, the Swiss mathematician 
Daniel Bernulli. A comprehensive theory of these sequen- 
ces is due to the great XVIII-th century mathematician, 
member of the Petersburg Academy academician Leonard 
Euler who devoted to the recursion sequences (series) 
Chapter 13 of his /ntroduction to the Analysis of Infini- 
tesimals (1748). Among the recent works the courses of 
the calculus of finite differences read by the renowned 
Russian mathematicians academicians P.L. Chebyshev and 
A. A. Markov should be especially noted; the courses 
include notes on the theory of recursion sequences. 


“ For the reader who is versed in analysis it can be stated that 
the theory is a strict analogue of the theory of linear differential 
equations with constant coefficients. 
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1. The concept of a recursion sequence is a broad gene- 
ralization of the concept of an arithmetic or geometric 
progression. Its particular cases also include the sequences 
of squares or cubes of natural numbers, the sequences of 
digits of a rational repeating decimal (and any recurrent 
sequence in general), the sequences of coefficients in a 
quotient of two polynomials in ascending powers of y etc. 
It follows that in the course of mathematical studies at 
a secondary school recursion sequences are very often met 
with. The theory of recursion sequences is a special part 
of the so called calculus of finite differences. This theory 
will be expounded here for readers with no previous 
experience in the subject (in only one place we shall 
refer without proof to a general proposition established 
by the theory of linear algebraic equations). 


2. Sequences will be written in the form 
Ws. Uy bag: aes Ue, Genes (1) 


or, shortly, {u,}. If there is a natural number ek, and 
numbers a,, a,, ..., a, (real or imaginary), such that 
starting from a certain number n and for all following 
numbers 


UntpH=UUny prt AUllnrp-ot--- +a,u, (Namo>1), (2) 


the sequence (1) is called a recursion sequence of order k 
and the relation (2), a recursion relation of order k. 

Thus, a recursion sequence is characterized by each of 
its terms (starting from a certain one) being expressed by 
the same number k of terms immediately preceding it, 
according to formula (2). The very name “recursion se- 
quence” is used just because one must turn back to the pre- 
ceding terms in order to obtain the following one. Some 
examples of recursion sequences will be given below. 

Example 1. Geometric Progression. For the geometric 
progression 


U,=a, U,=—ag, Uus=ag?, ..., u,=agr, ...3; (3) 
relation (2) takes the form 
Un+y = Uy (4) 


where k=1 and a,=g. Thus the geometric progression is 
a recursion sequence of the first order. 
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Example 2. Arithmetic Progression. For the arithmetic 
progression 


U,==a, Us=at+d, u;=a-+2d, ..., u,=at(n—l)d, 


we have 
Unyy =U, +d. 


This relation is not of the form of relation (2). Characte- 
ristic of the latter is that its right-hand member contains 
only sequence terms with constant coefficients. However, 
considering two relations for two consecutive values of n: 


Unse=Ungr td and u,y.= by, d 


and subtracting the latter from the former term by term 
we obtain 
Unte—Un41 = Un41 4s 
or 
Un+e2 = 2Uy41— Uy. (d) 


Relation (5) is of the form of (2). This suggests that the 
arithmetic progression is a recursion sequence of the second 
order. 

Example 3. Consider the ancient problem of Fibonacci * 
about the number of rabbits. Required is the number of 
pairs of adult rabbits descending from one pair in the 
course of a year if it is known that each adult pair of 
rabbits gives birth to a new pair every month, the new- 
born becoming fully adult in a month. The interest of 
the problem lies not in its solution, which is easily arrived 
at, but in the sequence whose terms express the total 
number of adult pairs of rabbits at the initial moment 
(u,), in a month (u,), in two months (u,) and in general 
after n months (u,,,). Obviously u,=1. One newborn pair 
will be added a month later, but the number of adult 
pairs will be the same: u,=1. In two months (from the 
initial moment) the little rabbits will become adult and 
the total number of adult pairs will equal 2: u,=2. Assume 
that we have calculated the number of adult pairs after 


* Fibonacci, or Leonardo da Pisa, a medieval Italian mathematician 
(circa 1170), wrote the book Liber Abact, containing extensive arith- 
metic and geometric information borrowed from peoples of Middle 
Asia and Byzantium and creatively reworked and developed by him. 


8 


r—I1 months, u,, and after n months, u,,,. By this time 
the u, adult pairs that existed before will have an off- 
spring of uw, additional pairs, so that after n--1 months 
the total number of adult pairs will be 


Unte=Ungy tUy. (6) 
Hence 
U,=u,+u,=3, UU oe O, 
u,=u,+u,=8, u,=U,--u, = 13, eee 


Thus the sequence obtained is 
Wey tg hy Ue 2 te, 
W205 Me By. US. 24 ey 


(7) 
i.e. every following term is equal to the sum of two pre- 
ceding terms. This sequence is known as the Fibonacci 
sequence (series) and its terms, as Fibonacci numbers. 
Relation (6) shows that the Fibonacci sequence is a recur- 
sion sequence of the second order. 

Example 4. Consider now the sequence of squares of 
natural numbers: 


eS A A, bee Uae 5 ae (8) 
where u,.,=(n+1)?=n?+2n+1. Thus 
Ung = u,+2, +1. (9) 
Increasing n by unity we have 
oe Upa ga nad: (10) 


And hence [subtracting (9) from (10) term by term] 


Unsgbng y= Unsy— Uy +2, 


or 
Une = 2Uy4y— Unt 2. (11) 
Increasing in equality (11) n by unity gives 
Un+3 = Uy yg—Unyi t 2. (12) 


Hence [subtracting (11) from (12) term by term] 
Un+g—Unse = LU yy — SU gay tun, 
or 
Un +3 = SU 4 2—SUnys Ug. (13) 
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The recursion relation obtained is of. the third order. 
Hence sequence (8) is a recursion sequence of the third 
order. Similarly it can be verified that the sequence of 
the cubes of natural numbers 


17, 23, 3%) ..., n3, .n. (14) 


is a recursion sequence of the fourth order. Its terms 
satisfy the equation 


Unga = 4g 3— OU ggg bt 4g yg 1 — Uns (15) 


and the reader is invited to derive it. 

Example 5. All recurrent sequences are recursion sequences. 
Consider for instance the sequence of digits of the rational 
repeating decimal 

761 


aang 0,57132132132 ... 


In this sequence 
U0; Uy 1, Ue, _— 


16 
He 2. ie 1 gy, She (16) 


It is obvious that 
Unyg HU, (n > 3). (17) 

Rewriting equation (17) to give it the form of relation 

(2) we have 
Ung =OUnge tO-Ugys +1 -U,. 

Thus, Eq. (17) is a recursion relation of the third order 
(R=3, a,=0, a,=0, a;=—1), and sequence (16) is a recur- 
sion sequence of the third order. 

Example 6. Now let us consider the sequence whose 
terms are the coefficients of the quotient of two polyno- 
mials arranged in ascending powers of x. Given 

P(x)=A,+A,x+...+A,x! 
and 
Q(x)=B,+Bx+... + B,x* (By 0). 

Let us divide P(x) by Q(x); if there is a remainder 
the division can be infinitely continued. The quotient will 
contain terms in the following order: 

Do+D,x*+D,x7?+ D3x°+...+D,x"+... 
10 


Consider the sequence 
=D Ue Dis. sig =D. ees (18) 


We shall prove that this is a sequence of order k (remem- 
ber that k is the power of the divisor). In order to do it 
we choose an arbitrary natural number n satisfying only 
one condition n>/—k+1, and in dividing we stop at 
the term of the quotient with x"+*. The remainder will 
be a polynomial R(x) containing powers of x higher than 
n-+k. Writing down the relation between dividend, divisor, 
quotient and remainder we obtain the following identity: 
Ajgt...+A,xt= 

= (Bot... + Byx*) (Dot... + Dua ax"**) +R (x). 


Find the coefficients of x**+* in the right-hand and left- 
hand members of the identity and equate them. As 
n+k>1+1, the coefficient of x"*+* in the left-hand member 
is zero. Hence the coefficient of «"** in the right-hand 
member is also zero. But the terms with x"** are present 
only in the product (B,+...+8B,x*)-(D,+...+D,4,*"**). 
(The remainder R(x), as was mentioned above, contains 
only higher powers of x). Thus the required coefficient is 

DiseBo tDrop-18i, +... +D,B ,- (19) 


According to the foregoing argument this coefficient is 
equal to zero: 


DisrBot+Dase-1Bit ar ari + D,,B,=9. 


Hence (as B,=« 0) 
B B 
Daik=— B Dnek-1— ++» —B Dn (1 >Il—k-+1). (20) 


This is a recursion relation of order k. It follows that 
sequence (18) is a recursion sequence of order R. 


3. Of all the above examples Example 6 is of the most 
general character. We shall now show that any arbitrary 
recursion sequence of the k-th order 


Una re tnt ioe. anes (21) 

satisfying the equation of the form 
Unte = AUnsper t+» FA Up (1>my> l), (22) 
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coincides with the sequence whose terms are the coeffi- 
cients of the quotient of a polynomial P(x) by the poly- 
nomial 


Q (x)= 1—a,x—...—a,x*. (23) 
Let n be an arbitrary natural number satisfying the 
condition n >rk+m—2; multiplying the polynomial Q (x) 
by u,tu,x+u,?+...+u,4,x%" we obtain 
(1—a,x—a,x?— ...—a,x*) (u, +u,x+ 
bee tpt tt tnt b 2. bye X") = 
= [uy + (tp—ayt,) HF. + 


Gg aes (Ups m—1 — UUs m—-2— aes —A,Um—1) eee ar 
ai (peep axa os —AyUm) xktm~*) al 
Hoses + (yg Og + — OU ge psi) XJ — 


—— [(Qiun a4 TF Ses + Aplln—p+2) xed ies + Agley 4X" 7"). (24) 


The polynomial in the first pair of square brackets is 
of a degree not higher than 1=k-+m—2; the coefficients 
in it are independent of the chosen number n. Lei us 
call it P(x): 


P (x) =u, + (u,—a@,u,) x+ 
Be pas eo ha a 20) 


All the coefficients in the polynomial in the next pair 
of brackets are zero by Eq. (22). Lastly, the third pair 
of brackets embraces a polynomial whose coefficients de- 
pend on n; it contains no terms of a power lower than 
n+l. 
Denoting this polynomial by R,(x) we can rewrite 
identity (24) in the following form: 


P (x) = (1—a,x—a,x?— ...— a,x") (u, + ux + 
fee Eiger t) +R, (2). (26) 
Hence, it is obvious that u,+u,x+...+u,,,x" is the 
quotient and R,, (x), the remainder of the division of P (x) by 
Q (x) = 1—a,x—a,x? —...—a,x* 
i. e. 
Wig Dee be say by Disa: ahd 


is indeed the sequence of coefficients of the quotient of 
polynomial (25) by (23). 
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Consider, for example, the Fibonacci sequence 
Ba. We i 2, Uy, UO oak 
Since its terms satisfy the relation 
Unse=Unsi tun (n > 1), 


then m=1, kR=2, a,=1, a,=1 and Q(x)=1—x—x’. 
The highest power in P(x) must not exceed k+-m— 
—2=1. By formula (25) we have 


P(x)=14-(—1Ix])x=1. 
Thus, the Fibonacci numbers coincide with the terms 


of the sequence of coefficients of the quotient of 1 by 
1—x— x’, 


4. One of the problems which are offered for solution 
to the students of secondary schools when studying the 
arithmetic and geometric progressions and the sequence of 
squares of natural numbers is to find the sum of n terms 
of each of these sequences. 

Let 


Up be: ecieige Ws nes (27) 


be a recursion sequence of order k whose terms satisfy 
the relation 
Unte=Qunrp-1 + AgUy 4 p—-9+ se + AU, (n = m). (28) 


Consider a new sequence formed by the sums s, of 
numbers (27); 


S, = Uy, S,==U,+ug, ..., S,=uUtugt... tu, ... (29) 


We shall show that the sequence of these sums is also 
a recursion sequence of order k-+1 and its terms satisfy 
the relation 


Sn+pti= (1 + a,) Snok 
+ (Ag—4)) Sn peat ++ F(A —Og-1) Sut —2p8n- (80) 
Note that 
eS Si. hg Sy HSS a. ss 
3] 
bay U,=S,—(Uy ae. +Uy_1)=S,—Sy—1) ea ( ) 
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Assuming s,=0 so that u,=s,—s, and substituting in 
equation (28) for u,, u,, ..., u, their expressions in terms 
Of Ss), S,, .--, S, We obtain 


Sntk—Sn+b—1 = (Sn4p—-1—Sntp—2) + 
+ a, (Sn4p—a7—Sn4e—3) + aia as + ap (Sp—S,-1); 


whence 


Sn+h= (1 +41) Sua p—1 + (Ga—41) Sng gat ++ + 
+ (A, —4p-1) Sp —4pSn-1 (1 Sm), 


or replacing n by n+ 1 we obtain 


Sntkea = (1+ 4,) Sy4p t+ (4g —GQ) Sngg-a te. + 
+ (Qg—Ay~1) Sn 41 — Ag Sp (n = m—1). 


This is a recursion relation of order k-++1. Several 
examples are given below. 
(a) Geometric Progression. In a geometric progression 


u,=—ag"-* and s,=u,+u,+...+u,=atag+...ag"™’. 
Since the terms of {u,} satisfy the relation u,,,=4qu,, 
the terms of {s,} must satisfy the relation 


Sn42= (1+ 9) Sn+17 9Sn- (32) 
(b) The Sequence of Squares of Natural Numbers. Here 

u,=n? and s,=—1+2?+...-+n?. Since the terms of {u,,} 
satisfy the relation 

unis = 3Un 42 —SUn 41 +U, 
(cf. p. 9) 

Sn +4 = 4Sn45— O8y49 +45n41—Sn- 
(c) Fibonacci Numbers. Since the Fibonacci numbers satisfy 
the relation 
Unte=Unsy tun 


their sums s, must satisfy the relation 
Sig = 25745 Si 


5. In the case of the simplest recursion sequences, for 
instance, arithmetic and geometric progressions, sequences 
of cubes or squares of natural numbers, and recursion se- 
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quences, any term can be found without calculating the 
preceding terms. In the case of the Fibonacci number se- 
quence or of the general sequence of coefficients in the 
quotient of two polynomials it seems, at first sight, to 
-be impossible to do that. In order to calculate, say, the 
thirteenth Fibonacci number we preliminarily find all 
the preceding terms (using the relation u,,,=U,+,+4,): 


Wa 1;. UH 1». Ug 2, 8. >», US: 
u,=13, us=21, uy= 34, u,,=55, u,, = 89, 
ies 144, te = 233: 


We now turn to a detailed study of the structure of 
recursion sequence terms. This study results in formulas 
which make it possible to calculate in the most general 
case any terms of a recursion sequence without resorting 
to the calculation of the preceding terms. These formulas 
can be regarded as far reaching generalizations of the for- 
mulas expressing the general term of an arithmetic or geo- 
metric progression. 

Consider the following recursion relation of order k: 


Unt e= AUy tpi Uglasp—aT +++ Tagua. (33) 
If it holds for any value of natural numbers n=1, 2, 
3, ..., then assuming n=1 we can obtain 
Ups, =AUp + AU,_ 1+... +Q,Uy,. 
Hence, knowing uw,, u,, ..., U,, we can now find u,,,. 


Assuming n=2 in relation (33) we obtain 
Ung, = AU yyy PAU, +... + Aglle. 


So now we know the value of u,,, too. In general if m 
is a natural number and we have calculated the terms 


Uy, Usy vvey Upy Upgyiy «++ Unep-1 


of a sequence, then assuming n=m in relation (33) we 
find the next term u,,,,. 

Thus the terms of a recursion sequence of order & which 
satisfy relation (33) are unambiguously determined by this 
equation if the first k terms of the sequence u,, u,, ..., U, 
are known. Choosing the sequences in various ways (there 
are no conditions limiting this choice) we can obtain an 
infinite number of different sequences which satisfy rela- 


15 


tion (33). The distinction between them will manifest 
itself in the first k terms (at least in one of them) and 
will be observed in the following terms too. 

For example, the relation of the first order 


Ung 1 = Guy, 


is satisfied by any of the geometric progressions whose 
common ratio is g (they differ in the value of the first 
term u,); the relation of the second order 


Unite = 2Un41—Un 


(OF Uy, 4o—Uny1 =Un4,—U,) is satisfied by any of the arith- 
metic progressions which differ in at least one of the terms 
u,=a and u,=a+d, i.e. differ either in the value of 
the first term (a) or in the value of the common diffe- 
rence (d) or in both. 

Now consider another relation of the second order 


Unieg= U, +1 lbs 


Besides the Fibonacci sequence 
diy: PegioG xs ae Oy: Oy LO, Dy OF. let 


where u,=u,=1 this equation is satisfied by many other 
sequences obtained by sampling various values of uw, and u,. 
For example, with u,=—3 and u,=1 we obtain the fol- 
lowing sequence: 


—3, 1, —2, —l, —3, —4, —7, —I1l, —18, —29, ... 


Let us have a certain number of sequences which sa- 
tisfy the same relation (33) 


Nig Mae, ead ands 
Ue Vas arty Udy saws (34) 
Z1y Lar sey ens oe @#y 


The following equations will hold: 
Xnth=UAntR-1 FAX nips as ee TP A,Xn, 
YntkeFUYntk-1 + QYn+p-2a ee + AY ns (35) 


enth = Uensp-i a ea pee ee + Bp2,- 


Taking as many arbitrary numbers A, B, ..., C, as 
sequences (34), we multiply each term of the first of equa- 
tions (35) by A, of the second, by B, ..., of the last, 
by C and add the products: 

AXnyn+ BYnsat e + C2nsn = 
=, (AXnspeat BYnge-at ++ + + C2n 44-1) + 
+ a, (AXn4%-— By nse =P ours eet oe, ae 
Gi aks + a, (Ax, a By, a are +Cz,). (36) 


It follows that the sequence 
t} = Ax,+By,+... +C2,, 


t,=Ax,+By,+...4+Cz,, 
Sch: de. de ok (37) 


obtained from sequences (34) by multiplying each term of 
the first sequence by A, of the second by B, ..., of the 
last by C and then adding the products termwise (first 
terms with the first ones, second terms with the second 
ones etc.) satisfies relation (33). Since the choice of the 
numbers A, B ...,C is arbitrary we can by changing 
them obtain, generally speaking, different values of the 
VOL MNS 05 Lag bas ae ave 
Now let the sequence 


Dg Ml sai ME oh wl fai (38) 
be such that it satisfies relation (33); the question is: can 
we give the numbers A, B, ..., C the values such that 


the first k terms of sequence (37) coincide with the first k 
terms of sequence (38)? If this can be achieved then accord- 
ing to the preceding argument all the terms of sequences 
(37) and (38) will also coincide, i.e. we will have for 
any natural number n: 


u, = Ax, +By,+...+Cz, (39) 


Thus it becomes possible (as yet hypothetically) to ex- 
press any of the infinite number of sequences which satisfy 
the same recursion relation of order k& in terms of some of 
sequences (34), by formula (39). The realization of this 
possibility depends on whether we can select the numbers 
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A, B,...,C so that the equations 
Ax,+ By,+...+Cz,=4,, 


Migs BUST: Ne aa: (40) 


Ax, 4 By, 7 £04 Up 


with arbitrary given values of their right-hand members u,, 
U,, ..., Uy, be satisfied. 

Since the numbers A, B, ...,C are in this case the un- 
knowns and the number of equations is equal to the k-th 
order of the recursion relation, it follows, that the number 
of the unknowns A, B,...,C [which coincides with the 
number of sequences (34)] should also be taken equal to &. 
It is known that a system of & algebraic equations (40) 
in Rk unknowns A, B, , C can have solutions depending 
on the kind of the coefficients in the system: X,, Y,, -- 

ie Cie wovaiy Dope: Ups cuy egy lo Gem the: kind ot the first 
terms of eaenes (34). The solutions will certainly exist, 
with arbitrary right-hand members uw,, u,, ..., u, if we 
assuine, for instance, that 


Reals Ye, 2iay 20; 


Ae w=! er ae (41) 


=o ti: ree 


Indeed, in this case system (40) takes a very simple 
form, and the solution is directly obvious: 


A=U,, 
DSU, 
Cu; 


Of course a different choice of numbers is possible with 
which system (40) has a solution whatever the right-hand 
members of the equations. Assume, for instance, 


Ki), Bet 5 ney eyed 
eat eae , 24; (42) 
en) =) » Zp= 1 
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The system now takes the form 
A+B+...+C=u,, 
B+...+C=u,, | 
C2 U;. 
Hence we obtain 
C=u,, ...,. B=u,—u,, A=u,—U,. 


We now turn to the general case and state the following 
theorem: 

For the system of k linear algebraic equations (40) in k 
unknowns to have one and only one solution A, B, ..., C 
with any values of the right-hand members u,, u,, ..., U, 
it is necessary and sufficient that the corresponding homo- 
geneous system 


Ax,+ By,+...+Cz,=0, 
ne iea oct e 40" 
Ax, + By,+ ...+Cz,=0 
have one only zero solution: * 
ASB = ys. 6 0 


The reader will easily ascertain that the condition in 
the formulation of the theorem is fulfilled in the particular 


* This is a convenient proposition since its application does not 
require any knowledge of the theory of determinants. The reader who 
is familiar with the theory will remember that for system (40) to 
have a solution with any values of the right-hand members of these 
equations it is necessary and sufficient that the determinant of the 
system 


My yw. By 
A Xo Yo 22 
a) Pa T 


be different from zero. The same condition is necessary and sufficient 
for system (40) to have one and only one solution for any given 
right-hand members (say, equal to zero). Thus for a system of & linear 
equations in Rk unknowns the condition for the existence of a solution 
with any right-hand members coincides with the condition that one 
and only one solution with zero right-hand members exists. Just this 
fact is expressed by the statement in the text. 
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cases (41) and (42). Later on we shall meet with cases in 
which the proposition stated above will prove useful. Until 
then we shall simply make use of the fact (established 
Se of the enw that there alwavs exist numbers 


,X,,.-., 2, [the first terms of sequence (34)] 
fe that ‘the. system of eigen (40) has a solution with 
any values of u,, u,, ..., U,- 


If such numbers have been chosen as the first terms of 
sequences (34), then according to the foregoing any sequence 
which satisfies the recursion relation (33) can be expressed 
by formula (39) in which the numbers A, B, ..., C are 
found from equations (40). The system of k sequences (34) 
by means of which the terms of any sequence satisfying 
relation (33) can be expressed by formulas (39) (i. e. by 
multiplying the terms by certain numbers A, B, ..., C, 
and adding the products) is called the basis of the recur- 
sion sequence. 

It follows from the above sequence that every recursion 
sequence has a basis which can be selected in different ways. 

For example, systems whose first terms are 


1 Oe. earns AD aes mre | 
0, le ’ 0 0, I, ’ l 
, OF 
05,0; aa 0;. 0; , 1 
(X) (k) 


form the basis of an arbitrary recursion sequence of the 
k-th order. 

Let us sum up the results of 5. 

For every recursion relation of order k there exists an 
infinit’ number of different sequences which satisfy this 
relation. Any such sequence can be formed from k sequen- 
ces which satisfy the relation and form its basis by multip- 
lying each of the k-th sequences by certain corresponding 
numbers A, B, ...,C respectively and adding the products 
ter mwise. 

Therefore in order to obtain a full solution of a recur- 
sion relation of order k it is sufficient to find but a finite 
number k of the sequences that satisfy it and form tts basis. 

We now give some examples in order to make this for- 
mulation clear. 
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Example 1. Let we have the following recursion relation 
of the second order: 


Use = 2Uy 91 — Ug. 


Its basis must consist of two sequences: 
xis Xo, X3, ee ey Xn vey 
Yi, Yo, Y3; ee ey Gn) siete 


In selecting them we shall assume 
Lal, ve) and 4. 0.<9,— 1. 


Since rewriting the recursion relation in the form 
Untea— 4n41 = Unt1—4ns 


we see that the difference of two successive terms of the 
sequence is constant, i. e. the sequence which satisfies this 
relation, is necessarily an arithmetic progression, then the 
sequence {x,' whose first terms are x,=1 and x,=1 is an 
arithmetic progression: with zero common difference, i. e. 


ide cee tees Heat) 


and the sequence {y,} with y,=0 and y,=1 is an arith- 
metic progression with the common difference equal to 
unity, i.e. 
0,152; co Aelene. (Y= n—1), 
A term of any recursion sequence which satisfies this 
equation can be expressed in the form 
u,= Ax, + By,=A + B(n—1) 
where A and B can be determined from the equations 
u,=A+B(1—)), 
u,= A+ B(2—1), 
or 
i= A, 
u,=A+B. 
Hence 
A=u,, B=u,—u, 
it follows that 
Uy, =U, + (n—1) (u,—ty). 
21 


This is the general formula for a term of any recursion 
sequence satisfying the relation 


Ue 2 pe ia 
Denoting u,—a, u,—u,==d we can give it the following 
form: 

u,=at(n—l)d. 


This is the well known formula for the general term of 
the arithmetic progression. 

Example 2. Consider another recursion relation of the 
second order: 


Unre are Uns =a U,. 


Assuming x,=1, x,=1 we obtain the familiar Fibonacci 
sequence: 
TOO Oy aa 


As a second sequence entering into the basis let us take 
the sequence {y,} in which y,=0 and y,=1. We can write 
then 


Ys=YotH=!1, W=HHty=2, ¥e=Wty,=3,... 


In this sequence y,=X,, Yg=%X_, Y,=X3, Yg =X, ..- and in 
general y, =X,-,(n=2, 3, ...). Indeed, since we have found 
that these equalities hold for any value of n<cm+1 so 
that in particular Y,,,,—=Xmn) Ym=Xm-1, then for y,,,. we 
shall have 


Yinse=Ym4i1 +Ym=XntXn-1 = Xin + 


and the equalities also hold for n=m-+2. 
Thus 


Yeats (2.05 cee): 
Therefore for any sequence, which satisfies the relation 
Unge = Ung1 Uys 
we find [formula (39)}]: 
u,, = Ax, +By,, 
where A and B are found from the equations: 
u, = Ax,+ By, =A, 
u,=Ax,+By,=A+B. 
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Hence 
A=u,, B=u,—u, 
and 
U, = U,X, ci (U,— Uy) Yn: 
With n>2 x,,_, can be substituted for y,; then 
Uy = UX, + (Ug— Uy) X,-4 (n > 2); 
or 


u,=u, (X,—X, -1) =F UX _-1° 
With n>3 


Xy =H Xya-r TX y-9s Or Xn Xn =) = Xn -25 
and consequently 
U,, = UX, -2 + UnX, —-1 (n = 3). 


Thus terms of any sequence {u,$, which satisfy the re- 
lation 
Unse = Uns F Uns 
can be expressed in terms of Fibonacci numbers by the 


formula we have established. In particular if u,=—83, 
u,=1 (cf. p. 16) then 


U,= — 3X,,~ FX y-1 (n > 3). 
6. It will be shown now that under some very general 
conditions the basis of the recursion relation (33) 
Unte = AUnsp-1 + Aguas p—2 zs orp a Apu ns 


comprising k geometric progressions with different common 
ratios can be found. In order to do it let us determine 
under what conditions the geometric progression 


X=1, %=G,...,%,= 9g", ... (GA) 
satisfies ralation (33). Noting that 
Xntk = a aa XntRe-1 ges, ieee xn = ges, 


we substitute these quantities in equation (33) (for u,4,, 
Unep—i» +++, &,) and obtain: 


grre-? = 09°" 2 ag" te 3 ae Aa g?s* 
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whence 
g*=a,g*-*>+a,g?-*+ ...+4,. (43) 


Thus a geometric progression can satisfy the recursion rela- 
tion (33) of order k if and only if the common ratio q of 
the progression satisfies the algebraic equation (43) of de- 
gree k with the same coefficients as those in equation (83). 
Equation (43) is called the characteristic equation for 
the recursion relation (33). If g=a is a root of the cha- 
racteristic equation (real or imaginary) then assuming 


Rear") Aneel Ze ime); (44) 


we obtain a geometric progression whose first term is 
X,==1.and with a common ratio a; this progression satis- 
fies relation (33). Since @ is a root of equation (43) we 
have, indeed, 


C= 0,0 9 062 2 OU, 


Multiplying the two members of this expression by a7! 
where n is an arbitrary natural number we obtam which 


Ce =a os OS ee sami ee eat +a,a"*, 


means that sequence (44) satisfies relation (33). 

Thus to each root g=a of the characteristic equation (43) 
there is a corresponding geometric progression (44) witha 
common ratio a; this progression satisfies the recursion 
relation (33). 

To form the basis of geometric progressions with 
different common ratios only one must have a sufficient 
number k of them and consequently have k different roots 
of the characteristic equation. 

Assume that all the roots of the characteristic equation 
are different: 


Gy=%, Go=B, ..-,9,=Y, 


We have then k geometric progressions which satisfy 
relation (33): 


Vis. Oe OL? eee Ole 5 ening 

2 n-1l 
oe (45) 
I, Vs sae ae 
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We shall show now that the system of sequences (45) is 
the basis of relation (33), i.e. that for any sequence {u,,} 
which satisfies relation (33) we can find numbers A, B,..., 
such that for any n we shall have 


u, = Aa"-'+ BBrs+ 2. + Cyt, (46) 
In order to prove this it suffices to verify that the sys- 
tem of equations obtained from (45) with n=1, 2,...,R 
A+B4+...4+C=u,, " 
Aa + BBA + Cy ty 


(47) 


has a solution with ines to the ne As Be hse ls 
with any values of the right-hand members of these equa- 
tions; for this it is sufficient, in turn, (cf. p. 19) that 
the corresponding homogeneous system 
A+B+...+C=0, | 
At BB .+Cy=0, (48) 


Acct! 4. BBE... Cpt =0 


have one and only one zero solution. That is actually so. 

In fact, suppose there exists a zero solution of (48), i.e 
there are numbers A, B, ..., C such that at least one of 
them, say, A is different from zero, and the numbers sa- 
tisfy system (48). To reduce this proposition to a contra- 
diction let us first form a polynomial M (x) of degree kR—1 


which vanishes with x=, ...,x=y and equals unity 
with x=a. Since this polynomial is of degree kR—1 and 
vanishes with R—1 different values of x: B, ..., y, it must 


be of the following form: 
M (x) =p (x—B) ... («—4), 


where p is some number. Taking x=-a, we must obtain 
M (a) =1; hence 


=p (a—B)... (a—y) 
and 
Y= (a—p)...(a—y)" 
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Thus 
__ («—B) ... («—Y) . 
a ~ (a—B)...(a—y)? 


and obviously the polynomial satisfies the given conditions. 
Removing the parentheses and collecting like terms we can 
give the polynomial the following form: 


M (x)=m,+mx+...+m,_,xF-?, 


If we multiply now equations (48) by m,, m,, ..., m,-, 
and add them term by term we obtain 


A(my+ma+... + m,_,a*-*)+ 

+ B (my + mB +... + my_,B4-") +... + 

+C (my +imyy +... +my_yv-4) =0, 
or 

AM (a)+ BM (B)+ ...+CM (y) =0. 
Since M(a)=1, M(B)=0,..., M(y)=0 it follows that 

A=0, 
which contradicts our assumption. 
Thus system (48) has one and only one zero solution 

and system (47) must have a solution (one only) with any 
U,, Uz, ..., u,, and this, in turn, means, that system (45) 


is the basis of recursion relation (83). 
Thus we have found that for every recursion relation 


Unt p= WUnrp-1 =F eee + a,U,, 
whose corresponding characteristic equation 


g* =a,g*-) +a,g'-? +... +4, 


has different roots: g=a, g=f, ...,g=y, there exists a 
basis consisting of & geometric progressions with the com- 
mon ratios: a, B, ..., y. In other words for terms of any 
sequence {u,,} which satisfies relation (33) there are k num- 
bers: A, B,...,C [they are found from equations (47)| 
such that 


= AGI BBY tae Cyr? (n=l, 2,0, ¢e2)s 
Let us sum up the ideas of 6. 
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A recursion relation of order k has a corresponding alge- 
braic equation of degree k with the same coefficients — its 
characteristic equation. Each root of the characteristic equa- 
tion is the common ratio of a geometric progression which 
satisfies the recursion relation. If all roots of the charac- 
teristic equation are different, R geometric progressions are 
obtained which constitute the basis of the recursion relation. 
It follows that in this case terms of any sequence which 
satisfies the recursion relation can be obtained by adding 
term by term certain geometric progressions (k in number ). 


7. We now turn to the application of the results ob- 
tained. Consider first the Fibonacci sequence. Its recursion 
relation is as follows: | 


Unse= Unay a, 


and the characteristic equation (43) accordingly takes the 
form 


g=q+l. 
Solving this equation we obtain two different real roots: 


lL lars a 


The general term of the Fibonacci sequence can there- 
fore be written in the following form: 
u,= Aa"-14 BBRrm?, 


In order to find the unknown coefficients A and B 
assume n=1 and n=2; we obtain 


u,=1=A+B 


uy=1=Aa-+ Bp=+(A + By +45 (A—B). 


Solving this system we find: 
om V5+1 B= V5—1 


2Y5 ' 9V5 | 
Hence 
y= VEE (14 VEY EVE (tye 
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Or 
(SEY), 


This is precisely the general expression for the Fibonacci 
numbers. At first sight the formula obtained looks unwieldy 
and poorly suited for calculations. However using it one 
can obtain many interesting results. We shall show, for 
instance, that the sum of squares of two adjacent Fibonacci 
numbers is a Fibonacci. number as well. 

Indeed 


a-$[(42)"+ (2-2 
tv E [CE (SYEY 2} 


Hence 


Une + Ua = 


— & 2 


It follows that 


Wnsa + Un =Uonsr (50) 
For example 


Uy, =U? + ui = 13? + 8? = 233. 


By the way, this is the answer to the Fibonacci problem. 
We leave it for the reader to prove a more general rela- 
tion for the Fibonacci numbers than (50), viz. 


Unlla tUnr Unt =Unemer- (51) 


The following theorem will be proved as another appli- 
cation of formula (49): 

For any two natural numbers a and b witha<), the 
number of consecutive processes in the Euclidean algorithm 
which are necessary to find the greatest common divisor 
(G.C.D.) of a and b is not higher than five times the num- 
ber of digits in a written down using the decimal system. 
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Applying the Euclidean algorithm to find the G.C.D. 
of the numbers b and a we obtain a chain of equalities: 


(1) b=ax'+y’, } 
(2) a=y'x"+y’, | 

3 = ie ad eae 

(3) yo=yX™ ry | (52) 
(B) yfBo2 = yD), | 

(RAL) yt) = yx (ht, J 


The consecutive remainders satisfy the inequalities 
Ay Sy Sy Se: Sy Dn > yy? >1. 

The remainder in the last of equalities (52) is zero. 

This means that the foregoing remainder y‘* is the re- 
quired G.C.D. of the numbers b and a. Therefore Rk is the 
number of operations required for finding the G.C.D. Our 
problem, as we have said, consists in evaluating the num- 
ber k. For this purpose let us compare the numbers y, 
y*-), ..., y’, a with the Fibonacci numbers u,, u,, us... . 
Note that y” >1=u, but the foregoing remainder y“-» 
is greater than y), and hence y“*-» >2=u,. Therefore 
we conclude from the equality & that 


2! yr + y* = ye 1), | + y) = is + U. — U, 
Thus 
y™ = is, yf i> ye 2) = Uy. 


Assume that we have proved the following inequalities 
to be valid 


GO SU ge YO Sipe UO Sas 
(m—1 > 2). 
Then from the equality y’”?~® = y@""-x" + y™ we obtain 
Yn ye Se Ye peel peage = Up a 
Continuing this argument we come to the inequalities 
Yo Zug, Yo Ups 
and further from (52, 2) we deduce that 
a=yx'+yY Sy AF Dugyr tug =Unse- 
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By formula (49) the term u,,, takes the form 


mpl FEY] 


Therefore - - 
sol (tE)—-(S78)"I> 
yee) 
(since od <1, and consequently Aas. oe 1). 
From (53) we find that 
(ste <aVb+1<V(a+1) < (YS) ayy 
(V5< Gt es a , as Vi< 3). 


Therefore 


€2345 eee: ' (54) 


We now note that 


an bade( (LE) -(I5¥)]'< 


Hence 


Consequently 
‘e \5k 
10k < 245 <(a+1). (55) 
If number a is written down in the decimal system with n 
digits (a is an n-digit number) it is obvious that 
10"7-2<a < 10"; 


hence 
a+1<10" 


and thus by virtue of inequality (55) 
10? < (a+ 1)®< 10" 
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or 
R< 5a. (06 


This is precisely the required result: the number of con- 
secutive k divisions in the Euclidean algorithm is less than 
five times the number of digits in the smaller of the two 
numbers, 6, and a, written in the decimal system. The 
foregoing proof shows that the most unfavourable case of 
the application of the Euclidean algorithm (as to the num- 
ber of operations near the limit established by the theorem) 
is that in which 6 and a are the adjacent Fibonacci num- 
bers. In order to confirm this statement let us take, for 
instance, D=u,,=6765 and a=u,,=4181. Here a is a 
four-digit nurnber and by the above theorem the number 
of operations in the Euclidean algorithm must be less than 
oxX4=20. In fact we have to perform in this case R=17 
operations: 


(1) 6765=4181x%1+2584 (10) 89=55x1434 
(2) 4181 = 2584%1+1597 (11) 55=34x1421 
(3) 2584=1597x1-+ 987 (12) 34=21x14+13 
(4) 1597= 987x1+ 610 (13) 21=13x1+ 8 
(5) 987= 610x1+4+ 377 (14) 13= 8x14 5 
(6) 610= 377x1+ 233 (15) 8= 5x14 3 
(7) 377= 233x1+ 144 (16) 5= 3x14 2 
(8) 233= 144x1+ 89 (17) 3= 2x14 1 


(9) 144= 89x1l+ 55 (18) 2 


Fibonacci numbers in descending order are obtained here 
one after another as the remainders of division. All quoti- 
ents (except the last) are unity and this is why the number 
of operations is so great. The greatest common divisor 
proved to be unity [equality (17)] which could be foreseen 
beforehand for the adjacent Fibonacci numbers. Indeed, 
from u,,,—=U,,,+u, it follows that the G.C.D. of the 
numbers u,,, and u,,, is the same as the G.C.D. of u,,, 
and u,. Therefore for any pair of adjacent Fibonacci num- 
bers the G.C.D. is the same. In order to find it one can 
just consider the pair u,=u,=1 which shows that the 
G.C.D. is unity. 


8. Let us consider as the next example the recurrent 

sequence (16): 
U0, Ue 1s ely. UO oy. Up, U5. aes 
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The recursion relation in this case is 
Unygeu, (n > 3), 


and thus the characteristic equation is 


g = 1. 

The roots of this equation are 
_ — 1,,;V3 ee wee 
=I. p=—5+t-5— and ‘eae eee a mae 


Therefore the general term of the sequence should be as- 
sumed to have the following form [cf. formula (46)]: 


u, = Aa"-14 BRr-14 Cyn-t — 
=A+B(—44i — ‘+¢(— s—! ae 


We can require this formula to be valid for all values 
of n with which our recursion relation holds too: n=3, 
4, 9, 


Note that 
| V3 ALF oe ee 
aa se a — (cos $—i sin} . 
l - V3 ae Bi» tat, TE 
ti a [cos 5 +i sin 5 | 


Therefore by De Moivre’s theorem (formula) 
uy, =A+B(—44i aS “+e(— si a = 
=A+(—1)""1B [cos 3 (n— ) isin fo) 
+(—1)*-?C jcos $ (n—1) +isin F (n—1)] = 
=A+(B+C)(—1)""! cos 3 (n—1) + 
+i(—B+C)(—1)""! sin} (n—1). 


Let us have B+C=A, and i(—B+C)=A,; then the 
above formula will be rewritten in the form 


u, = A+ A, (—1)""! cos = (n—1) + A,(—1)"™ sin} (n—1) 
(n> 3), 
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and we have now to find only the unknown coefficients 
A, A, and A,. Taking n=3, n=4 and n=09 we obtain 
three equations in three unknowns: 


; 3 

u,=1=A-+A, cos +. A, sin = A—> A, + — ee Ax, 
u,=3=A—A, cos — A, sn 7 =A+A,, 

u,=2=A +A, cos + A, sin 242 Ae, 


Hence we find 


A=2, A,=1 and Aaa = 


V3 ° 
Thus 
u,z=2+(—1)"™ hay (n—l) Fe sin (n —1) >| — 
= 2+(—1)". ie sin(n—2)4 (n>3). 


It can be seen that in this example the general term of 
the sequence is expressed by trigonometric functions which 
is in agreement with the recurrent character of the sequence. 

As the last example consider the one directly related to 
the division of polynomials. 

Given two polynomials: P(x)=3-+x«?--x®5 and Q(x) = 
==2—x—2x?-+ x3. The problem is to find the structure of 
the coefficients in the quotient of P(x) by Q(x). The se- 
quence of the coefficients in the quotient 


= Di USD ig. ei lg a 


as noted in 2, is a recursion sequence whose terms satisfy 
equation (20): 


B B 
Dase= ay Di epsi ah caer D,, (n SIl—R+ I), 


where k& is the degree of Q(x), B,, B,, ..., B, are the 
coefficients of Q(x) and / is the degree of P(x). 
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This means that kR=3, B,=2, B,=—1, B, =—2, B,=1, 
{=65: 


2 l 
Di+2 +5 Dasi—z Dn (n Z=5—3+1 = 3), 


et 


l 
Diis=F Di +a + Dnti—z Pn (n > 3). 


The characteristic equation takes the form 
a I 
T= 77 1I—-F: 
or 
3 l 2 pene l \ —_— 
P—9— 3 (F—1) = (9—z) G—N Gt) =0 
Therefore its roots are 


] 
a= 3 B=1, yvy=—l, 


and for D, we obtain the formula 
D,= A(z x) +B +C(-1" (n>3), 


Assuming n=3, n=4, and n=5 we obtain three equa- 
tions: 


D,=4A+B-—C, 
De=GA+BHC, 
D,=yA+B—C. 


In these equations not only the coefficients A, B, and C, 
but also the numbers D,, D,, D, are unknown. In order to 
find the unknowns let us perform the division of P(x) by 
Q (x) up to the point where the quotient will contain terms 
of the fifth power, inclusive. 
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Dividing we get 
3+ x2 — x5 

_ 3 r 
3— > X— 3x? + 5x8 


| 2—x—2x?2+ x3 
3 5 3 5 3 19 
SRLS Wed Fa ec Vc eh 


3 3 
ane 2. O¥y8__ +5 
5 + 4x 9 * x 


3 3 3 3 
Seah | 2 3 4 
x x ae ar 


2 4 
42 e—S ts 
45 x29 x84 x8 +25 x8 
25 x84 4x4 3= x8 
“oe 15 xt 28 xe 4 1 x! 
Bax lame 
© BS xt 95 85S x 4 De x’ 
Lo x8} Act O28 27, 
Hence 
D,=2, D=4, D,=25, D,=15, D,=2, =H 
and the system of equations obtained above takes the form 
ZA+B—CH15, 
KA+BLC=25, 
+ A+B—C=2, 


whence we find 


Thus 
Dit oe at GaS3 
n=4e-satataglcl)* jes). 
The problem is solved. It follows by the above formula 


51 179 179 
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9. In all of the above examples the characteristic equa- 
tion had only simple roots. Now consider the example 
pertaining to the sequence whose terms are the sums of 
squares of natural numbers (given on p. 15). The recursion 
relation for this sequence is 


Sn+a = ASn +3 — OSn 42+ 48y41—Sns 


and accordingly its characteristic equation takes the form 


q* = 49° — 6q? + 4q—1, 
or 
q* —4q° + 69° —49 + 1 =(q—1)* = 0. 


It has only one quadruple root: g=1; therefore we ob- 
tain in this case only one geometric progression with the 
common ratio 1 whose terms satisfy the above recursion 
relation. 

In such cases we have to find other simple recursion 
sequences which can form, together with this geometric 
progression, the basis of our equation. Such sequences in 
our example are: 


(the reader can easily verify this assertion). Without treating 

this very general case which would require rather unwieldy 

calculations we shall consider the following typical example. 
Let the recursion relation be 


ney CE Ottty p-y— Ch Oty pat ».. +(—1)P Charu, 
(57) 
where C&-1, Ck-?, ..., C2 are binomial coefficients of the 


k-th order. The corresponding characteristic equation 
gt = Ch tagt4—Chtaagt*# +... + (1) Chat 
can be written in the form 
(q—a)*¥ = 0. 
It has a k-fold root g=a; obviously 
(a—a)* = ak —Ch ak + Ch*at— ... +(—1)*Chak=0. (58) 
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Consider in general the following identities 
(a—a)P-™ — oF —m Ci i hak 


+ Choe tote (1) "Ch_,a*" = 0, 
where m=O, 1, 2, ..., R—I1 or 
(1—1)*- —m __ = Crm m he baa —M—2_ es 


+ (—1)PCRom (1) Ch_, = 0. (59) 
Equality (59) corresponding to m=0 takes the form 
Ch—CEI 4 CEA... + (1h CY + 


4... 4(—1)CR=0 er? 
Noting that 
Pm 32 (ns) ere (Ves bd ( ee Re (cee oe ae 8 ee ee 
Sia ee (kR—p)— (R—m—p+l)... (R—p) Chom 
(m=1, 2,...,kR—1; OS pak—m), 


OF 
k(k—1)... (R—m-+1) Cit = 
=(k—m—p+1)...(k—p) Ci, (60) 


we multiply each of the equalities (59’) by the correspon- 
ding multiplier k(R—1)...(k—m-+1). Wecan now using (60) 
write them in the form 


(k—m+1)... RCE—(k—m) ... (R—-1) CEO+ 0+ 
+ (—1)"(k—m—p+ 1)... (k—p) Ce” + 
+ o.. $(—1)™ 1-2... mCh=0 = (59") 
(m=1, 2,..., R—1). 


We shall now prove that the following equations hold 
with m=0, l, 2, ..., R—I: 
h™Ch—(k—1)™ CEA... (1) # (Rp) CRY 
+...+(—1)*-07Cf=0. (61) 
Indeed, the equality corresponding to m=0 coincide 
with (59’) and thus holds true. 
Assume, reasoning by induction, that equalities (61) have 
already been proved for m=O, 1, ..., j ({<Rk—2). Let 


us demonstrate that the equality which corresponds to 
m=j-+1 also holds true. In order to prove it we shall 
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make use of a polynomial of degree j+ 1: 


f(x) =(@—f) (@—j41)...@—-1) x= 


=x Byxi—... BX. (62) 
Multiplying (61) with m=1, 2, ..., j by the correspond- 
ing numbers 6, B,, ..., B; we obtain 
B,kCE—B, (R—1) Ch + ) 


+... +81 (—1)* (kp) Ce? + 
+...+8, (—1)*-0-CZ=0, 
gn eh de why eS, te. BB (63) 
BkICE—B, (k—1 VCE + 
$0 $B (1)! (hwy CE + 
4 $B, (—1)#-0-C2=0. J 
Let us rewrite equality (59”) corresponding to the value 
m=j+1 in the form 
f (k) ChR—f (R—1) Ce*+ ... +(—1)* f (Rp) CE + 
+...+(—1)*f(0)CR=0. (64) 
In doing that we made use of the equalities 
(k=). R= F(A), 
(k—j—l)...(ek—l)=f(k—1), ..., (R—p—jJ)... 
...(R—j)=f (R—j),.... 
Adding (63) and (64) term by term we have: 
[Bikt ...+B,A/+f (k)] Ce— 
—[B.(2—1) +... +B, (&—1 +f (R—1)] C7 +... + 
+ (—1)* [B, (Rp) +... +B; (Ap) +f (A—p)] CR + 
foe. (—1)* BOF... +B,-OF+ f(0)] CR=0 
But by (62) we can write 
Bix t+Bix?+ 2... +B pel + f (x) = xt}. 
Therefore the above result yields 
k/**CE—(R—1)/** CE + 2. + (—1)" (R—p)itt CR + 
foe + (1-0/2. CR = 0, 
This is precisely equality (61) for m=j+1. Thus we 
have demonstrated that Eqs. (61) hold true. 
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Lastly, consider an arbitrary polynomial whose degree 
is not higher than k—1: 


Py Ag Aga Fe a eA. (65) 


Multiplying equalities (61) with m=0, 1, 2,..., R—1 
by corresponding factors A,, A,, ..., A,_, we obtain: 


A,Ci—A,Ce* + oa i 
+ (—1)" A,CeY+ ... +(—1)* A,CZ=0, 
A, RCE—A,(R—1)CE'+...+ 
+(—1)" A, (k— 1) CE H+... +(—1)* A,-0-C2=0, 
Ag_sk®'CE—Ag_, (R—1)8 7 CE +... + 
4+(—1)" A,_, (R—pyt CEH Ee + 
hfe Ag, OF CRO, 


Adding the equalities term by term we obtain 


(A, +A,k+...+A,_,Rk*-') CE— 
—[A, +A, (R—1) +... + Agia (R—1)*-*] Ch + 6 + 
“ar ot Ay (Rp) +... + Ag-s (R—-py3] ChB + 
$e. (—DAA, + Ar 0+... + Aga OF] Ch = 0, 


Or 


P (k)-Ck—P (k—1)-Ck2+4....-+(—1)# P(0)-C2=0. (66) 


And hence an arbitrary polynomial P(x) whose degree 
is not higher than R—l1 satisfies relation (66). 

Assume, in particular, P(x)=(x+n—1)", where n is 
an arbitrary natural number and m, an integer OSmM< 
<k—l. 

Then equality (66) takes the form 
(k+-n—1)"™*Ck—(k+n—2)*CPI+...4 

+(—I)*(n—1)" Ch=0, 


or multiplying it by a**”-? and substituting unity for C%: 
(k+n—1)" aktn-1 — Ch 0 (k+n—2)™ htt 2 
— Cot (kR4+n—3)"oktn-84 0 4 

++ (—1)F-* Chat (n—1)™ a"-?. (67) 
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Comparing (67) with (57) we see that each of the fol- 
lowing Rk sequences satisfies the recursion relation (57): 


Le Oy Os. oe ay OOS nas (m= 0); \ 
O, a, 2a%,..., (n—l1) a}, ..., (m= 1); | 
0, a, 22a?,..., (n—1)?at?,..., i= 2); (68) 


Oy3O. BP ar ne Ce Pra a. 4.45% Setna east, | 
If we come to establish that these equalities form a ba- 


sis it will follow that the general term of any sequence 
which satisfies equation (57) has the form 


Uy = [Bo +B, (n—1) +... + Byoy (n— 1) 4] a"? = 
=Q(n—1)a"-4, (69) 
where Q (x)=B,+B8,x+...+8,_,x*71 is a polynomial of 
a degree not higher than k—1 with arbitrary coefficients. 
In order to prove that (68) is a basis it is sufficient to 
demonstrate that the following system of k linear equations 
Bv,+8,:0+...+8,_,:0% t=u,, 
B,+B8,-:14+...+8,_,: 17 !=4,, 


B,+B,(k—1)+...+B,_,(k—1)*"! =u, 


has a solution with respect to the unknowns B,, B,, ..., Bs_; 
with any u,,..., 4, (which follows from the proposition 
stated on p. 19), i.e. that the system 


B,=0, 
By,t+8,+ o oe +B,_,=0, 


© ee @e@ @ @ @ @  @ @ @ @  @ @ @  @  @ @  @  @ 


Baba 1) 8 ot dO Be = 0 


has only a zero solution. But the equations of the latter 
system show that 


Q(0)=Q (I=... =Q(k—-1)=0, 
i.e. the equation 
Bvt Bix-+...+B,x%"-*>=0 
of a degree not higher than k—1 has at least & different 
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roots: 0, 1, 2,..., R—1. Hence 
BD =p, Sn. = Bea 0, 


and thus we have proved that sequences (68) form the 
basis of the recursion relations which satisfy equation (57). 
If the recursion sequence satisfies the general equation 


Uns p= Uns per $F AUasp—2t +» + Aglly (4,79), (70) 
its characteristic equation 


g®=a,gP-14+...-+a, (71) 
can have a root a of multiplicity a, a root B of multip- 
licity b, ..., a root y of multiplicity c, and a total of 


at+b+...+c=k roots. 
In this most general case it can be demonstrated that 
the basis is formed by the following k sequences: 
| RY 3 Fk! © “nc © Aaa ee ee 
O20, gen (ISO 4, 25 
5» PS Oak cane [2 yo eos 
p 


Oye he (Ns ts 


Therefore 

u, = Q(n—1)a®-?+ R(n—1)B"-*+ ... + S(n—1)y""*, (72) 
where Q(x), R(x), ..., S(x) are some fixed polynomials 
of a degree not higher thana—1, b—1, ..., c—1, respec- 
tively. 


It follows that the general term u, of any recursion 
sequence has the form of a sum of the products of polyno- 
tials with respect to (n—1) (or, which is the same with 
respect to n) by general terms of geometric progressions 
whose common ratios are equal to the roots of the characte- 
ristic equation (71). 

If all the roots of the latter equation are simple, the 
above polynomials are constants and the general term of 


4] 


the recursion sequence takes the form of a sum of the terms 
of geometric progressions. 

It can be proved that the converse also holds true, na- 
mely, that any sequence {u,} whose general term is express- 
ed by formula (72) is a recursion sequence *. 

The corresponding equation (71) is constructed using its 
roots a, B, ..., y and their multiplicities a, b, ...,c (which 
are the increased by unity degrees of the polynomials Q, 
R, ...,S). Hence the recursion relation (Eq. 70) can be 
found immediately. 

Consider as an example the sequence 


u, = (n—1)2-27-2 4 37-2 


Comparing it with (72) we see that the roots of the 
characteristic equation are: a=2, B=3, a@ being a mul- 
tiple root of multiplicity 2+ 1=3. Therefore the charac- 
teristic equation must take the following form: 


(q—2)? (q—3) = q* —9q? + 30q?— 44q + 24=0, 
and the recursion relation is 
Un +4 = IUn43— 30U, 4. + 44u,,,—24u 


mn’ 


It is left for the reader to ascertain that the above se- 
quence satisfies the last equation. 


10. We shall illustrate now the results of section 9 by 
several examples. It was shown in section 2 that the terms 
of an arithmetic progression satisfy the relation 


Untg = 2Uyn41—Un, 
the squares of natural numbers, the relation 


Untg= Stings — Sug a1 aS Uns 


and the cubes of natural numbers, the relation 
Unsa= 4€Uyp3— Sg ge +4 41 — Uy. 


Obviously all these equations are particular cases of the 
equation 


* The proof of the theorem and its converse, which we omit, can 
be found in Chapter 4 of the book Division which Leaves a Remainder 
in Arithmetic and Algebra, Moscow-Leningrad, 1949; the theory of 
recursion sequences is expounded in this book in a manner different 
from that used in this booklet. 
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Una = CR Uy a p-1— Ce Uy gpa be» (TPO Chu, (87') 


discussed in section 9 (in this case a=1). 
The general term of any sequence which satisfies this 
equation must have the following form [cf. formula (69)]: 


u, = B,+B,(n—1)-+...+B,_,(n—1)*?-*. — (69") 


In order to find the coefficients B,, B,, ... B,_, it is 
sufficient to solve the following system of & linear algebraic 
equations in & unknowns: 

By =U), 


B,t+8,+...+8,-,=,, (73) 


By +B, (k—I) + +2. + Byoy (R— 1)? = wy. 


In the case of an arithmetic progression k=2 formula 
(69’) takes the form 


u,=B,+-B, (n—1)), 
and system (73), the form 
By=u,, B,+B,=4,. 


It can be seen that B,=u, is the first term of the pro- 
gression and B, =u,—u,=d, its common difference. It fol- 
lows that 

U, =u, +d (n—l) 


which is a well known formula. 

There is no need in making calculations for the case of 
sequences of squares or cubes of natural numbers since we 
know beforehand that the general terms are u, =n? or 
u,—=n*. However, there is some interest in applying the 
relations (69’) and (73) to the derivation of formulas for 
the sum of an arithmetic progression, as well as for the 
sum of squares or cubes of natural numbers. 

In section 4 we proved that if the terms of a sequence 
{u,} satisfy the equation of the form 


Unth=UWUnsp-1 tT Wu4nsre-2 ++ Tagua, 
then the sums {s,,$ of the terms of these sequences (s,=4,, 
S,=U,+U,, S;=U,+Uu,+U,, ...) Satisfy the equation 
Sntker = (1 +4) Sang + (GQ: —Q) Suspea bese + 
+ (Ay— p= 1) Sn 41 — 2p Sp- 
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Obviously for equation (57’) we have 
a,=Ch, a,=—Ch, ..., a,=(—1)*-*CE. 


Therefore 
I+a,=1+C,=Chi, 
Pe cee h+ Ch) + — Char, 
ea ee 


—A,= (—1)*#Ck= (ech 


and the equation for {s,} can be written in the following 
form 
Snenti = ChiSn+e—ChuSntp-at--- + (—1)*C#fis,,, 
or 
Sn+eti— CheiSn+e + ChsiSn+e-1— = + (—1)*** Chis, = 0. 


Thus, if a sequence {u,} satisfies the recursion relation 
(57') of the k-th order, then the sequence of the corres- 
ponding sums {s,} satisfies a relation of the same form 
but of order k+1. In particular, for the arithmetic pro- 
gression k=2, for the sequence of squares of natural num- 
bers R=8 and for the sequence of their cubes R=4 and 
it follows that for the sequences of the respective sums we 
must take k in the above equalities (57’), (69’), (73) great- 
er by unity: 3, 4 and 5. 

(a) The Sum of the Terms of an Arithmetic Progression. 
According to what has been said s, is expressed by for- 
mula (69’) (substituting s, for u,) with k=3. It follows 
that 


= B,+ B,; (n—1)+ B, (n—*” 


The coefficients B,, B,, B,, are dete rom system 
(73) (also substituting s, for uw, with k- 


By =s, =U, 
B,+8,+8B8,=s,=u, +4, =2u, +d, 
B,+2-B8,+2°B,=s,=u,+u,--u,=3u,+4+ 3d. 


Solving this system of equations we obtain 
l l 
Be U.. B,=u,+ 54, B,= => 4. 
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Whence it follows that 
I I 
S,= Hy + (4 +34) (n—1)+34(n—lP = 
=nu, + 5d (n—1)n = Pet 


Mls eee +(n—1) dm (ts +n) 
2 2 


(b) The Sum of Squares of Natural Numbers. 
Taking in formulas (69’) and (73) R=4 and _ substitut- 
ing s, for u, we have 


S,=By+B, (n—1)+B, (n—1)?+B, (n —1)° 
and 
Bias, =. )5 
B,-+8,+ B,+B,;=s,=1--2?=5, 
B,+2B,+4B8,+8B,=s,=142?+3=14, 
B,+3B,+9B,-+27B,=s,=1+2?+3?+ 42=30, 


The last system yields 
By=1, B,=2-, By=15, By=z. 


Therefore 
13 3 : 1 5 
s,=1+ 3 (a—1)+5 (2-1) += (n—1) = 
n+son+—ni= ness en) an nine ene) . 


= 6 

We have obtained a well known formula. 

(c) The Sum of Cubes of Natural Numbers. This sum is 
expressed by the formula 

5, ME 

We leave the derivation of this formula for the reader 
as an exercise. 

In conclusion let us consider one more example, the 
sequence a, 207, 3a3, ..., na”, ... (240, al). 
In this sequence 

iene”. (a 1 28 inv) 
It is easily seen that 
Ung = 20g 4 — WU. 
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Indeed 
20Un, 41 —O7U, = 2a(n+ 1) a"*t?—ana” = 
= (n+ 2) Ort? — Un yo. 

Since R=2, a, =2a and a, =—a?, the sequence consist- 
ing of sums {s,} (s,=a, s, =a + 2a?, S,= a+ 202 + 
+ 3a’, ...) satisfies the equation [cf. (30)]: 

Sn43 = (a, +1) Sno2 t+ (@g—4,) Sagi — as, >= 
= (20+ 1) 5,4.—(a + 2a) Snti + 0?S,. 
The corresponding characteristic equation is as follows: 
g® = (20+ 1) g?— (a2 + 2a) q+ a 

It can be easily ascertained that this equation is satisfied 
with g=a. Dividing the polynomial g*—(2a+1)q? + 
+ (a?-+2a)q—a? by g—a we obtain the quotient 

qg?—(a+l)qt+a. 

For this reason the two remaining roots of the characte- 
ristic equation satisfy the following equation: 

g—(at+l)qt+a=0. 

These roots are @ and 1. 

Thus the characteristic equation has a multiple root a 
of multiplicity a=2 and a simple root B=1. 

Therefore we obtain for s, [cf. formula (69) in which we 
must write s, instead of u,, a=a, Q(x)=B,+B,x, i.e. 
a polynomial of degree 1, B=1 and R(x)=C,, a constant]: 

Ss, =(B, +B, (n—1))a"-'4+-C, (n=1, 2,3, ...). 

The coefficients B,, B, and C, are found from the fol- 
lowing system of equations with n=1, 2 and 3, respecti- 
vely: 

B,+C,=s,=%, (8, +8,)a+C,=s,=a+ 2a, 
(B,+2B,)a?+C, =s, =a + 2a? + 3a. 


Hence 
3__ 9g? 2 
B= Se, B= and Gao. 
hagas 


BP sb ecco a2 
(a—Ip 
— bn® 2 (Un 41—Un) 
(a — 1)? 
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CONCLUSION 


This book is meant to give the reader an idea of the 
variety of recursion sequences and of their place in mathe- 
matics. Along with it it was shown that the recursion 
sequences do not greatly differ from their simplest forms — 
the geometric progression and the sequences whose terms 
are powers of natural numbers (in particular, the sequences 
of natural numbers as such being an arithmetic progres- 
sion); the recursion sequences can be expressed making use 
of these simplest sequences. 

However even in elementary mathematics sequences are 
very often met with which are not recursion sequences. 
One of these is the sequence of prime numbers, one of the 
most important in mathematics: 


25950 Ty Lye AS, L7G). 235: exc 


This sequence with its deep and complex properties is 
a subject of the theory of numbers. 

Sequences of values of many elementary functions 
are mot recursion sequences as well. For instance, the 
sequence 

1 ft 4 l 
» 9? 3 4? ce ey n? cet 


l 
(this is the sequence of values of the function y=— with 
x=1, 2, 3, ...), or the sequences 


log!, log2, log3, log4, ..., logn, ... 


(these are the sequences of values of the functions x and 
log x) etc. 

These and similar sequences * (including the recursion 
sequences) are the subject of the branch of mathematics 
mentioned above—the calculus of finite differences. 


* Sequences of values of the so-called analytical functions are 
meant whose simplest instances are the elementary functions. 
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And finally, in elementary mathematics and particularly 
in analysis studied in colleges and universities a very 
important place is occupied by convergent sequences, i.e. 
sequences which have finite limits. Their study is the most 
important problem in the theory of limits and belongs to 
the foundations of analysis. The properties of individual 
terms of the sequences are even less than secondary: only 
the existence of a limit and its value are of consequence. 

We consider these remarks to be necessary in order to 
bring home to the reader that the theory of recursion 
sequences given in this book from the viewpoint of its 
subject-matter as well as from the viewpoint of the rela- 
tionships that result from it is only a particular and very 
modest chapter in the theory of sequences. 
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